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Abstract. The concept of effective temperatures in nonequilibrium systems is studied within 
an exactly solvable model of non-Markovian diffusion. The system is coupled to two heat baths 
which are kept at different temperatures: one ('fast') bath associated with an uncorrelated 
Gaussian noise and a second ('slow') bath with an exponential memory kernel. Various 
definitions of effective temperatures proposed in the literature are evaluated and compared. The 
range of validity of these definitions is discussed, ft is shown in particular, that the effective 
temperature defined from the fluctuation-dissipation relation mirrors the temperature of the 
slow bath in parameter regions corresponding to a separation of time scales. On the contrary, 
quasi-static and thermodynamic definitions of an effective temperature are found to display the 
temperature of the fast bath in most parameter regions. 



1. Introduction 

Recent years have witnessed enormous efforts in understanding glassy systems by studying a 
broad variety of model systems. At the same time, it has been found that various other systems 
like colloids, granular materials, etc. share several features of typical phenomenology of glassy 
systems. Therefore, a common macroscopic description of glassy systems would be extremely 
helpful. However, unlike equilibrium thermodynamics, there is no general, well-established 
framework for a macroscopic description available today that works out of equilibrium. 

A key role in the macroscopic description of glassy systems could be played by the concept of 
effective temperature Based on the fluctuation-dissipation relation (FDR) pQ, the effective 

temperature is supposed to govern the slow, not equilibrated degrees of freedom, while the short 
time dynamics and the equilibrated degrees of freedom are governed by the bath temperature. 
Subsequent molecular dynamics simulations |3J HJ on a supercooled model fluid during aging or 
in shear flow are in agreement with this picture. However, it has been argued that non-linear 
generalizations of the fluctuation-dissipation theorem (FDT) out of equilibrium could spoil the 
interpretation of effective temperatures from FDR 5 . The authors emphasize that "Physical 
interpretations of the several temperatures and their relation to experimental probes and to their 
mutual relation in non-equilibrium steady states are needed." (p. 2014 in [H]). In addition, more 
recent molecular dynamics simulations [6 j found that effective temperatures should sometimes be 
defined from static FDR rather than from the dynamical one, as was done in [HIE] • The authors 
of [0] conclude their study with a "puzzle: when should one use static linear response and when 
should one use a time-dependent relation?" . For a review on the numerical evidence of effective 
temperatures see also (Zj. Up to now, experimental investigations of effective temperatures in 



glassy systems are not conclusive. Experiments on spin glasses in the aging regime [H] and 
self-diffusion of granular material in a Couette cell [H] have shown some evidence in favor of 
the existence of an effective temperature as introduced in pQ. In dense colloidal suspensions, 
very recent experiments |lf)j found no deviation from the FDT , while an effective temperature 
could be measured in a colloidal glass of laponite I12j. Most experiments, however, involve 
frequencies that are rather high compared to the inverse relaxation time. Earlier experiments 
on an aging colloidal glass ^SI observed violations of the FDT in electrical but not in rheological 
studies. 

In view of these controversies, it seems appropriate to resort to simplified model systems, 
where the validity of effective temperature concepts can be studied in more detail. It has been 
found, that mean-field spin glasses with one-step replica symmetry breaking scheme do allow 
the definition of an effective temperature from the FDR |14j . while those with continuous replica 
symmetry breaking do not [7j. Field-theoretic calculations on the critical dynamics of spin 
models show, that the effective temperature defined from the FDR does depend on the variables 
considered |15l I16j . Also in the Ising chain with Glauber dynamics and dynamical trap models, 
a meaningful effective temperature cannot be defined from the FDR .17 . A slightly different 
perspective was taken in |181 1191 12U] . where Brownian particles coupled to two heat baths at 
different temperatures were studied. The aim of the present study is to derive and compare 
different static as well as dynamic definitions of effective temperatures in these systems and 
discuss their range of validity. 

This paper organized as follows. The model of non-Markovian diffusion in the presence 
of two heat baths is introduced in Sec. El In Sec. [3 the model is simplified by considering 
the overdamped limit. From the exact solution of the model, different definitions of effective 
temperatures are evaluated and compared in Sec. Some conclusions are offered in Sec. 03 

2. Model 

Consider a particle of mass m at position x with velocity v moving in a potential V(x) under 
the influence of two thermal baths. One bath is held at temperature T s i ow . Its influence on the 
dynamics of the particle is described by the retarded friction coefficient (memory kernel) T(t). 
The other bath is kept at temperature Tf ast . Contrary to the first, slow bath, the correlation time 
of the second bath is small enough, so that the particle experiences an instantaneous friction 
described by the friction coefficient To- The equations of motion read 

x = v 

dV f* 

mv = -—- J dsT(t-s)v(s)-T Q v(t)+Z(t)+ri(t). (1) 

The fast bath is modelled as Gaussian white noise with (r](t)) = and (rj(t)r](s)) = 2Tf ast To5(t — 
s), whereas the random force due to the slow bath is characterized by (£(i)£(.s)) = T s \ ow T(t — s). 
We choose units such that Boltzmann's constant = 1. The model (^Q) generalizes earlier work 
on non-Markovian diffusion [2J| by including a second heat bath at a different temperature. 

The diffusion equation (^Q) is difficult to study in general due to its non-Markovian character. 
For the special case of an exponential memory function, 

r(t) = -exp[-t/(a7)], (2) 
a 

it is possible to map the non-Markovian dynamics (^Q) onto a Markovian process in an extended 
state space [2T] . 



x = v 



dV 

mi = -— + z(t) -T v(t) + r](t) 

z = -l v (t)-—z(t) + C(t) (3) 
a cry 

where £ is a Gaussian white noise with 

{Cms)) = ^L S (t-s). (4) 

Equations © together with the noise correlators specify the model to be considered in the 
present study. 

3. The overdamped limit 

Formally, the model of non-Markovian diffusion introduced in Sec. |2] can be mapped onto a 
system of two coupled Brownian particles, each equipped with its own heat bath. In order to 
further simplify the analysis, we here consider the overdamped limit mv — * where the inertia 
term can be dropped. In the overdamped limit, the time evolution equations Q simplify to 

dV 

r x = -— -T(0)x(t)+T{t-t o )x(t )+rj(t) + h(t) 

h(t) = -J*ds dT{t g ~ S) x(s) + C(t) (5) 

These equations are studied also in [T^]. In ^S], the initial time is set to ~~ * — oo, thus 
T(t — to)x(to) — ► and the term depending on the initial condition is absent. 

Employing again the exponentially decaying memory function P|). closed, Markovian time 
evolution equations for x and h are obtained. In a generalized notation 1 , the resulting time 
evolution equations take the form 

1 dU 

xi = -=r"a h^lX 2 + 7?l(t) 

I i OXi 

X2 = -tt%2 + v 2 xi + r/ 2 (t) (6) 

1 2 

where the Gaussian white noise is characterized by (r]i(t)) = and (r]i(t)rjj(s)) = 2TiT~ 1 5ij8{t — 
s). In Eq. we have neglected the transient term T(t — to)x(to), which is justified for times 
t long enough or for the initial condition a: (to) = 0. Although Eqs. © are very similar to those 
studied in JH]) we point out the important difference that the memory function ((21) does not 
satisfy r(0 + ) = as is assumed in [T9] . 

Only for the special case T\V\/Ti = T2V2/T2 (in terms of the original parameters, this 
condition reduces to T s i ow = Tf ast ) is the equilibrium distribution function corresponding to 
© a Boltzmann distribution, p e q(xi,X2) oc exp {—TJ(x\)/T\ — Kix\jT2 + V1T1X1X2/T1). In 
that case, the marginal distribution f e q(xi) = f^° oc d,X2P e q(xi, x 2 ) is given by / eq (^i) oc 
exp (— U{xi)/Tx + K2\x\ /2T2) where k 2 \ = (r 2 ^2) 2 1&2- Thus, the coupling to x 2 leads to a 
shifted equilibrium distribution, that can be described by a repulsive harmonic potential of 
strength k 2 i. In case the potential U(x\) is itself harmonic, U{x\) = k,\x\/2, the shifted 
equilibrium distribution / eq is of the same form as the uncoupled equilibrium distribution but 
with a renormalized spring coefficient. If instead one insists on the equilibrium form /eq(^i) oc 

1 The original model is recovered if the following identifications are made: x\ = x, X2 = h, U = V + (2a)~ 1 x 2 , 
k 2 = a, I\ = r , T 2 = a 2 7, v\ = l/r , ^2 = l/(a 2 j), Ti = Tfast, T 2 = T s i ow . 



exp(— f7(xi)/T eq ), an effective temperature 2 is defined by T eq = T\/[\ — riT2^i ^2(^1^2) ]■ 
In terms of the original variables, U is harmonic if the potential V in Q is harmonic, 
V(x) = k\ i qx 2 /2. Thus, k± = k\ : q + q _1 and the corresponding identification of an effective 
temperature equals the bath temperature T eq = T\ as it should, since in this case the condition 
Ti^i/Ti = T2V2/T2 reduces to Tf ast = T s i ow . It is interesting to note, that the same effective 
temperature T eq is obtained by adiabatic elimination, if X2 is assumed to vary much faster than 
xi [2J. Then, the assumption K\/T\ <C K2/T2 replaces the condition Ti^/Ti = T2^2/T2- 

Yet another approach to obtain effective temperatures is to employ the Quasi-Equilibrium 
Approximation (QEA) which is used successfully in many areas of statistical physics 
125] . Following the standard procedure [12], extremizing the entropy functional S[p] = 
— Jdxidx2plnp subject to the constraints of fixed moments Mq = (1), M\ = (a? 2 ), M2 = 
(x 2 )> M3 = (X1X2), the quasi-equilibrium distribution is found to be given by p*{x\,X2) = 
exp (Ao + \\x\ + X2X2 + A3X1X2). The Lagrange multipliers Aj are chosen such that the 
constraints are satisfied identically. In equilibrium, p* reduces to p eq . Out of equilibrium, p* 
approximates the true non-equilibrium state by equilibrating all degrees of freedom except for the 
macroscopic variables Mj. The quasi-equilibrium entropy S* is defined by S*(M) = S[p*] and 
satisfies dS* = ^2,-XjdMj. Within the QEA, the temperature is defined as T^ A = dS*/dU. 
For simplicity, we will assume a harmonic potential, V(x) = ki$x 2 /2. Then, the effective 
temperature within the QEA is given by Tq^ a = (2/k\)Xi. Thanks to the simplicity of the 
model, the Lagrange multipliers can be evaluated explicitly by performing Gaussian integrals. 
The resulting expression for Tqea is 

Tqea = Ki((xl) - (xix 2 ) 2 I (x%)). (7) 

In the uncoupled case, (X1X2) = 0, Eq. (J7J) reduces to the equipartition theorem and Tqea = T\. 

In the following section, we evaluate and compare different definitions of effective 
temperatures derived from the exact solution in case of harmonic potentials. 

4. Exact solution for harmonic forces 

For the harmonic potential U(x\) = K\x\/2, k% = fc^o + l/a, Eqs. @ are stochastic differential 
equations of the narrow-sense linear type that can be solved exactly by transformation to normal 
coordinates. The solution reads 

Xi{t) = Yl KaX a (t) (8) 

where the coefficients bi >a are defined by = — a_/(a + — a_), b\- = a + /(a + — a_), 

62,+ = — o-)) &2,- = — — a -) and 

X a (t) = e^-^XM + / ' cfee-^-^TfcOO. (9) 

Jt 

The noise terms appearing in Eq. © are related to those in (JBJ by rj a (t) = rji(t) +a a rj2{t). The 
eigenvalues associated with Eqs. © are given by 

2 Note, that T eq ^ T\ if even both baths are at the same temperature T\ = T2 since the effective temperature 
mimics the effect of the coupling on the equilibrium state. Therefore, quantities that depend only on x\ but not 
on X2, show equilibrium expectation values corresponding to T oq rather than Ti. 



with D = (ki/Ti — K2/T2) 2 + kv\V2- The coefficients a± introduced above are defined by 
a± = (ki/Ti — c±)/i>2- Note, that the special case 1*2 = has to be treated separately. In 
that case, the time evolution equation for %2 can be solved independently of x\, resulting 
in an additional noise term in form of an Ornstein-Uhlenbeck process resulting from X2- In 
particular, one finds in this case, c± = K\/Ti and a± = —(ki/Ti — K2/T2) vi- In the uncoupled 
case, v\ = 0, the eigenvalues reduce to the individual relaxation frequencies c+ = K2/T2 and 
c°_ = K\/T\. In the coupled case, v\V2 7^ 0, stable solutions with positive eigenvalues c± exist for 
k-ik-2 > (ri^i)(r2i / 2)- For v\V2 < a further condition |«i/ri — K2/T2I > 1\J\v\Vi\ is necessary 
in order to ensure real values of c±. For the original parameters, v\ = 1/To > and the stability 
condition reads ma > 0. 

In the long-time ('static') limit, we find from the exact solution (JHJ) 

^ C a +C(3 \Li I 2/ 

In the uncoupled case, v\ = 0, the equipartition theorem h\ = T\ is recovered. Here and 

in the following, R\ = k\q refers to the bare spring constant of the original model. If, however, 
one considers Eqs. © as starting point, h% = k% is interpreted as the bare spring constant. 
Note, that in terms of the original model parameters, the uncoupled case v\ = corresponds 
to To, a — > 00. In the general case ^ / 0, one might assume that some generalization of the 
equipartition theorem holds, 

r static = K\ (xl)^ (12) 

with an effective ('static') temperature T stat i c . Note, however, that such an effective temperature 
results from an underlying canonical distribution function only for the special choice T\V\/Ti = 
Y2V2/T2, as was demonstrated in Sec. 01 For this choice of parameters, T sta ti c = T s \ ow . For 
a general choice of parameters, Tstatic 

denotes the effective temperature deduced from the 
equipartition theorem if we would not be aware of the coupling of x\ to X2- A comparison 
of these and other definitions of effective temperatures is provided later. 

Next, we consider possible definitions of effective temperatures from time-dependent FDR. 
The response of the system © to a time-dependent perturbation f(t) is measured by the 
response functions 



Ri,i(t, t) 



S( Xj {t)) f 



(13) 

/=o 



Sf(t' 

Due to the simplicity of the model, the response functions are also obtained analytically, 

R l , 1 (t,t') = T^ h^-^eit-t') (14) 

a=+- 

where 0(i) denotes the unit step function. The response functions are causal, R(t,t') = for 
t < t', and time-translational invariant, Ri t i(t,t') = Ri,i(t — t'), as they should. The correlation 
functions are defined by 

djit, = (xi(t)xjtf)) - (xi(t)) ( Xj (t')) (15) 
Inserting the exact solution, Eqs. (JHJ), ©, into Q15j). one obtains 

£7y (t, = V + aaa ^) L-^)*, _ e -(*-*)c e -(^)*l (16) 

c a + c p \T 1 T 2 J L J 



where t > t' has been assumed without loss of generality. Again, due to time-translational 
invariance Cij(t,t') = C^j{t — t') holds. For the special choice of parameters corresponding to 



the original model, Eqs. l[T%|) and (fl6|) are equivalent to the equations given in ^ 
the long time limit of the equal time correlation function reduces to lim^oo Cij(t, t) 



Note, that 

From Eq. (|16|) we observe that the decay of the correlation function is determined by the two 
eigenvalues c±. This can be interpreted as a two-step process, where the initial, fast decay is 
described by the larger eigenvalue c_, while the smaller eigenvalue c+ governs the long time 
decay. Such an interpretation is particularly relevant if the two eigenvalues are well-separated. 

The fluctuation-dissipation relation (FDR) is most conveniently expressed in terms of the 
integrated response function, defined by Xi,j{t) = ' ' 
differences one obtains 

-l \^ 



Xi,i(t) 



-II 



ds Ri j 

OO l 'J 



1 



(s). Integrating Eq. (|14|) over time 
a ) • (17) 



- tc, 



In equilibrium, one can prove under quite general assumptions the validity of the FDT , 

1 



Xi, j (t) = -[1-C ijj (t)] 



(18) 



where Xi,j(t) = Xi,j(t)/Cij(t = 0) and Cij(i) = Cij(t)/Cij(t = 0). In the present case, the 
FDT holds identically in the uncoupled case V\ = 0. For v\ ^ 0, however, the FDT is violated 
for the system studied here, as can be seen directly by comparing Eqs. (|T6|) and (|T7j). It was 
proposed in pQ, that a meaningful effective temperature can be defined from the FDR by 



- eS 



da 



(19) 



hi 



For short times, the effective temperature defined from Eq. ()19|) coincides with T%, as is readily 
shown by expanding the exponentials in Eqs. (|T6|) and (fTTjl to first order in tc a . For long times, 
the time dependence of both the correlation and response function is dominated by the smallest 
eigenvalue. In this regime, the effective temperature defined from Eq. (|19|) is given by 



plong 
-eff 



1 2 



«2 

r 2 



(20) 



Finally, we introduce also the effective temperature T a 



1 / c 



x (i ^ oo)/Ci,i(0), Ci,i(0) 



, introduced in [£]. This temperature, defined from the 'static' limit of the FDR , was found 
to be useful when the definition (|2()|) was problematic 0. From the long time limit of (|17j) we 
find 

Too = Static [l - r i r 2 VlV 2 (KiK 2 y 1 ] (21) 

where r sta tic is the static, effective temperature defined in ((T2*|) . Note, that is well-defined 
since we assumed k\k 2 > (Fi^i)(r2^2) above. Within the original model, = T sta tic holds 
which is directly verified inserting the original model parameters. 

Table ^ summarizes the different definitions of effective temperature considered here. 

In Fig.^ the integrated response Xi,i(*) is plotted in a parametric plot, versus the correlation 
function Ci ; i(t). Parameters are chosen as 7 = Tq = Tf ast = 1, Ki t o = T s \ ow = 2. The values of a 
are varied between 1 < a < 10. From Fig^one observes a cross-over from the short time regime 
with T e g = T\ to the long time regime with T e ff given by (|20|) . This scenario has been found also 
for spin glass systems ^ and in molecular dynamics simulations of sheared model-glasses 
For increasing coefficients a or /-t^o, T e g increases and the cross-over is shifted to later times. 
For increasing 7, T e fj increases as well, however, the cross-over point is not shifted. 



Table 1. Overview of different definitions of effective temperatures. 



symbol 


effective temperature from 


defined in 


Tqea 


quasi-equilibrium approximation 





^static 


equipartition theorem 


m 


T c s 


FDR 


Pi) 


rnlong 
1 cS 


long-time limit of the FDR 




T 

J- CO 


'static' limit of the FDR 





The time evolution of the effective temperature (|19|) is shown in Fig. [21 where T e g is calculated 
with the help of the time derivatives of Eqs. (|16ft and (|f Equilibrium initial conditions for the 
uncoupled system were chosen. Figure [2] shows the approach of the effective temperature to its 
long time asymptotic value (|19j) . For the present choice of parameters, the asymptotic value is 
approached for times t > 1, while T e fj remains close to the bath temperature T\ for short times. 

Figures |21 \M show a comparison of different definitions of effective temperatures for several 
model parameters. We observe that for increasing a, or 7 or Ki q, the static temperatures all 
approach the temperature Tf ast of the fast bath, while the effective temperature from the long- 
time asymptotics of the FDR approaches the temperature T s i ow of the slow bath. Only in case of 
very strong friction coefficient Tq associated with the fast bath, does T e g fail to approach T s i ow . 

5. Conclusion 

The concept of effective temperatures has been studied within an exactly solvable model of 
non-Markovian dynamics for a harmonic oscillator coupled to two heat baths held at different 
temperatures. For parameter ranges (a > 1, or 7 > f, or kiq 3> 1) corresponding to 
a separation of time scales of the fast and slow bath, it has been found that the effective 
temperature T^ ng defined from the long-time limit of the fluctuation-dissipation relation (FDR ) 
indeed approaches the temperature of the slow heat bath. This result is therefore consistent with 
the proposition [Q of T^ ng as a useful definition of effective temperatures out of equilibrium. 
For intermediate values of the parameters, when the separation of time scales is not achieved, 
the effective temperature T^ s is intermediate between the temperature of the slow and fast 
heat bath. Only in the unlikely case of a friction coefficient Tq associated with the fast bath, 
larger than the friction 7 associated with the slow bath, T^ 18 approaches the temperature of 
the fast bath, while the quasi-equilibrium temperature is somewhat higher. 

In comparison, other, static effective temperature definitions yield results that are much more 
difficult to rationalize. In most cases, the temperature of the fast bath is obtained from these 
definitions. However, in some cases, values intermediate between T s i ow and Tf ast are obtained 
even when time scales are well separated (see e.g. figures 3 and 4, panels b and c, for T stat i c or 
Tqea)- 

We finally mention that a breakdown of the effective temperature from the long-time limit 
of the FDR can be observed in the model (jHJ) for particular choices of the parameters that are, 
however, inconsistent with the original model |HI For couplings v\ < 0, v<i > 0, Y\ < T2 and 
Ti > T2, there are parameter ranges where an overshoot in the response function x 1S observed 
together with a corresponding undershoot in the correlation function C. There, the parametric 
plot x versus C continues to negative values of C and x(C) becomes multi-valued. Approaching 
this regime, the plot x versus C flattens more and more, which makes the application of 
(|I9[) difficult. A similar situation was observed for some observables in molecular dynamics 



simulations jSj. 

The present study confirms and extends previous results |23j that several, non-equivalent 
effective temperatures can be defined. In the model system considered here, which could be 
taken as a model for an internal degree of freedom in a slowly driven system, the definition 
that appears to have the most sensible behavior is associated with the long time limit of the 
fluctuation dissipation ratio. 
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C u /C u (0) C u /C u (0) 

Figure 1. The integrated response function ()17|) is plotted versus the correlation function (|16|) . 
Dashed lines have inverse slopes calculated from Eq. (|20|) . The inset shows the time evolution of 
the normalized correlation function for the same parameters. The parameters have been chosen 
as a = 7 = Tq = ki ; o = «a = 1, Tf ast = 1, and T s i ow = 2 if not stated otherwise. In the left 
figure, the values of a from bottom to top are a = 1, 2, 5, 10, while 7 varies from top to bottom 
on the far left of the right figure as 7 = 1,2, 10, 100. 




time 

Figure 2. The time evolution of the effective temperature defined from the time-dependent 
FDR, Eq. Q19JI . is shown on a logarithmic time scale. The same values of the parameters as in 
Fig. ^ (left panel) are used. Dashed lines are the long-time effective temperatures calculated 
from Eq. (|2"U|). 




Figure 4. Different definitions of effective temperatures are shown as a function of model 
parameters a, 7, To, and K\fi. The remaining parameters are chosen as in Fig. Q with a = 1 
and 7 = 100. 



